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TECHNICAL MEMORANDUM X-64907 


A STEERING LAW FOR A ROOF-TYPE CONFIGURATION FOR A 
SINGLE-GIMBAL CONTROL MOMENT GYRO SYSTEM 

INTRODUCTION 


As a candidate for the momentum exchange device to be used for attitude 
control of the Large Space Telescope (LST — a payload of the Space Shuttle) 
and the High Energy Astronomy Observatory (HEAO), four Single-Gimbal 
Control Moment Gyro (SGCMG) systems have been investigated recently by 
many researchers. 1,2,3 

Two configurations of four SGCMG systems have been proposed thus 
far: a pyramid-type configuration, i. e. , the gimbal axis is normal to the faces 
of a square- based pyramid ( Fig. l) , and a roof-type configuration ( Fig. 2) . 

At first the pyramid-type configuration was investigated extensively and 
various steering laws were proposed. In this configuration, however, the 
singular states (the state at which the torque output axes associated with each 
of the SGCMG' s in the system are coplanar and the system cannot respond to 
out-of-plane commands 2 ) cause serious difficulties and each of these proposed 
steering laws has shown some unsatisfactory performance in certain situations. 

As a configuration in which the difficulty of singularities would be less 
serious, the roof-type configuration was investigated and, at the same time, the 
OMEGA steering law was proposed. 2 The main feature of this configuration is 


1. B. G. Davis, A Comparison of CMG Steering Laws for High Energy Astron- 
omy Observatories (HEAO's), NASA TM X-64727, 1972. 

2. J. W. Grenshaw, 2-SPEED, A Single Gimble CMG Attitude Control System, 
TR- 243-1139, Northrop Services, Inc., 1972. In this reference this law is 
called 2-SPEED (Two Scissored Pair Ensemble, Explicit Distribution) , but it 
is usually referred to as OMEGA (Optimum Momentum Exchange by Gimbal 
Alignment) (see p. ii of this reference). 

3. J. R. Glease, A Summary Description of the 2-SPEED Steering Law and 
Configuration for Single Gimbal CMG's, NASA Internal Memo, S&E-AERO-DOI/ 
Analytical Investigations Section, October 1972. 




that it consists of two pairs of gyros, with each pair sharing a common gimbal 
axis direction. 

Although the OMEGA steering law gives rather good control performance, 
it still has the defect that stable singular states exist. A stable singular state 
is defined as a singular state which can be maintained by setting the input torque 
command to zero. The existence of such states is not desirable because of the 
following reason: If a certain torque command is given when the system is in, 
or near, a singular state, there will be large control deterioration, and the 
existence of stable singular states makes the possibility of occurrence of this 
control deterioration larger. 

In this report, for four SGCMG systems with the roof-type configuration, 
a steering law is proposed which does not have the above-mentioned defect. 

This steering law is obtained by regarding the CMG system as a sampled-data 
system and providing a new momentum distribution scheme. 

The basic procedure of the steering law is as follows: (l) at each sampl- 
ing instant the desired total momentum at the next instant is calculated from the 
present total momentum and torque command; ( 2) according the the predeter- 
mined momentum distribution scheme, the desired total momentum is distributed 
into two desired momenta for the two CMG pairs; (3) a desirable combination 
of gimbal angles which realizes these desired momenta is calculated for each 
CMG pair; (4) the desirable gimbal angles are compared with the present 
gimbal angles to select a best way to attain the desirable gimbal angles, and 
(5) a gimbal- rate command is calculated and if the calculated rate command 
exceeds the hardware limit, the rate command is modified not to exceed it. 

An important feature of this steering law is that, because of a new 
momentum distribution scheme, it does not result in any stable internal sin- 
gular state. Therefore, the possibility of the occurrence of unfavorable torque 
command when the CMG system is at, or near, a singular state is negligible. 
Moreover, it will be shown by computer simulations that even if such an unfor- 
tunate situation happened, control deterioration would be very small. Another 
feature is that, since this momentum distribution scheme treats the momentum 
directly, any reasonable momentum distribution can easily be realized. 

In the next section the roof-type SGCMG system is described and some 
considerations are given to the system from the viewpoint of sampled data 
control system. In the section titled Steering Law, a steering law is described 
in five steps. In the section, Momentum Distribution Scheme, such a scheme is 
given which is left undecided in the Steering Law section. The CMG-out operation 
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is discussed in the section having that title, and results of the computer simula- 
tions are given in the section, Computer Simulations. Redistribution out of 
undesirable initial states, insensitivity of the steering law to the sampling inter- 
val of the CMG system, the behavior of the CMG system after reaching any 
momentum saturation, and the flexibility of the steering law with respect to the 
desirable momentum distribution are considered in the Discussion section. The 
conclusions are given in the last section. 

In the following sections, the physical explanations and selection of con- 
stants for the simulations are made with regard to the LST. Almost the same 
argument can be presented for the HEAO. 


ROOF -TYPE CMG SYSTEM 


Figures 2 and 3 show the roof-type configuration of a four SGCMG 


system and its mounting arrangement 
relative to the LST. Each CMG angular 
momentum ( and torque) vector is 
restricted to the plane I or II, which is 
skewed relative to the vehicle Y-Z plane 
by the angle 0. . 

The angular momentum vector of 
the i-th CMG is denoted by h^ 1=1,2, 

3,4. In addition to the coordinate sys- 
tem (XYZ), we also use coordinate sys- 
tems (X x Y Zj) and (X n Y n Z R ) , which 

are shown in Figure 2, for the pairs (h , 


TELESCOPE AXIS 


X 



h ) and (h , h ) , respectively. 

ju O a 

The angle between the rotation 
axis of the i-th CMG and X for i = 1, 2 

( or X^ for i = 3,4) is called the i-th 

gimbal angle and is denoted by a: ^ (in 

degrees). Figure 2 shows the state 
where a. = 0, i = 1, 2, 3, 4. 



SOLAR WING 


Let us assume that h. ' s have the 
same magnitude h and define 


Figure 3. Mounting arrangement 
of roof-type CMG system 
relative to LST. 


z 


4 


Then, 


-I = th Il J h I2 | ° lT ’ ( inX i Y I Z I “ coordinate system) 

-II = th IIl ' h II2 ' 01X1 ( inX II Y II Z H " coordinate system) 
where the superscript T denotes transpose and 


h Il 

= h + (cos a + cos « 2 ) , 

(la) 

h I2 

= h^ ( sin a + sin a ) , 

(lb) 

h m 

= h + (cos a 3 '+ cos a 4 ) , 

(lc) 


and 

h n2 = h * ( sino; 3 + sina 4 ) * ( ld ) 


The total momentum h of the CMG system is given by 


4 

h = 2 h. = h + h 

— . „ —l —I —II 

i=l 


i 1 T 

= [h h , h ] , (in XYZ-coordinate system) 

x J y ; z 


( 2 ) 



where 


and 


h x = sin 0 ( h i2 + h E2 ) * 


h y = h n - h ni ’ 


h z = COS 0 (h I2 -h n2 ) ; 


h can also be expressed as 


h 



HI 


(in Yj. Y Y - coordinate system) . 


( 3 ) 


Notice that h is determined uniquely if a A [ j a 2 ! j a ^ ] is 

given, but that a is generally not determined uniquely even if h is given. 

The input (control) variable of the CMG system is the rate of the 
gimbal angle r (= a ) ; the output is the torque, i. e. , the rate of the total 
momentum t = t , and the purpose of this system is to make t_ equal to the 
torque command . 

For the case where a digital computer is used to obtain the desirable 
gimbal rate from the torque commandj,^ , it is impossible to make_t_ equal to 

t at all times. Hence, the treatment of the system as a sampled-data system 

— c 

is appropriate in this case. This treatment also plays an important role for the 
development of a steering law in later sections. Let t be the time variable, A 
be a sampling period, andh (t) and a (r) be _h and a at time r . 
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It will be simple and practical to impose the restriction, 

a (t) = a (nA ), nA = r <(n+l)A,n=0, 1, 2, ... 

This corresponds to using a zero-order hold. The output torque then satisfies 
the following equation: 

'j* ^ 

j t_ (r) dr = h [a (nA) + “ ( n <A) ( T * - n ^) 1 
nA 

- h [a (nA)] . ( 4 ) 


Since the torque command is usually obtained from the observed position 
and rate of the spacecraft and the spacecraft is exposed to unknown disturbances, 
it is impossible to know the future t (t) , t > nA exactly at the present time 

nA . Therefore, it is natural to assume that the purpose of the system is to 
make the average of the output torque during the coming A period, _t (r) , nA S 
r = (n + 1) A, equal to t (nA) where n = 0, 1, 2, . . . . 

The left-hand side of equation (4) is the momentum which is transferred 
from the CMG system to the spacecraft and which, in turn, causes the change 
of angular rate of the spacecraft. Hence, the above purpose of the system can 
also be interpreted as the effort required to change the rate of the spacecraft 
a specific amount during each sampling interval. 

Hereafter, a (nA) andh (nA) will be denoted by a (n) andh(n) for 
the sake of simplicity. 


STEERING LAW 


In this section a steering law is proposed which will fulfill the purpose 
described in the previous section. This steering law consists of the following 
five steps. Assume that the present time is r = nA . 

Step 1: Given the present gimbal angle a (n) and torque command 
t (n) , calculate the desired total momentum h , at the next sampling instant: 
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( 5 ) 


h^ = h (n) + A t^ (n) 


" [h dl | h d2 | h d3 J 


(in Yj Yjj Y-coordinate system) , 


where 


(in Yj Y^j Y-coordinate system) , 


h I1 (n) = [cosa 1 (n) + cos « 2 (n)] , (6a) 

h I2 (n) = h^ [ sin ( n) + sina 2 (n)] , (6b) 

h m (n) = [cosa 3 (n) + cos (n)] , (6c) 

and 

h JI2 (n) = h + [sinct 3 (n) + sina 4 (n)] . (6d) 


Step 2: Obtain the desired momentum distribution between two momen- 
tum vectors h^ and h^ at the next sampling interval: 

h Il = f I '-d ’ - <“)' ’ 


h (n) = 


h K (n) 


h n2 (n) 


h n (“) - h ni (n) 
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( 8 ) 


h ni = " f n [ V - (n)1 


h I2 h dl ’ 


and 


h II2 h d2 


( 9 ) 


( 10 ) 


f and fjj. are given by the momentum distribution scheme, which will be 
described in the next section. 

Remark 1: It is clear from the relation (3) that, for h^ and h^ to satisfy 


111 + -II — d * 


three relations, (9), (10) and 


h I2 h ni h d3 * 


( 11 ) 


should hold. Therefore, all that the momentum distribution scheme can do is 
to distribute h among h and h . This distribution is expressed as a pair 

UiJ ixj- 

of functions and f^. 

Step 3: Calculate the pairs of gimbal angles (“n’ a u> and ("o' a m ) 
which give the momentum vectors h^ and h^, respectively ( Fig. 4) . As can be 
seen easily from Figure 4, a ^ anda I2 are given by 


a 


II 


y i * 5 i 


(12a) 


and 


\ 


9 



(12b) 


O' 


12 



6 


1 


where 


>1 



- tan -1 (h n /h I2 ) + 90 sgn (h J2 ) 
90 sgn (h^) [ sgn (h^) - ll » 

90 - tan 1 [h I /(4h s|e 2 - h^ 2 )] 
0 , ifh t ^ 2h* , 


, ifh 12 * o 


ifh I2 = 0 > 


, ifhj < 2 h^ 


(13) 


(14) 


and sgn ( . ) is a sign function defined by 


sgn (a) 


+ 1 , if a > 0 

0 , if a = 0 

- 1 , if a < 0 . 


Similar equations give a ando^. 

Remark 2: The magnitude of h , h j should be less than or equal to 2 h* 

in order to be physically realizable by the pair I. However, there is a possibility 

that h may become larger than 2 h . The inequality sign in equation (14) takes 

I * 

this into consideration. 

Step 4: For pair I, select the better way between the following: (a ) to 
bring to and o 2 to or (bj) to bring to a J2 and ce g to a^. 
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2 2 2 2 

Then we decide as follows: Choose (a^ if e^ + 5 e I3 + e J4 ; choose 

(bj) otherwise. 

For pair n, define (a n ) , (b n ) , e m ~ e^, and follow the same pro- 
cedure. 

Remark 3. It may be better to choose (a^) if max ( t ©j I » I e gl ) = max 
( I e 3 l » I e 4 l ) and choose (bj) otherwise, but the above decision method was 
selected because of its simplicity. 

Step 5; Calculate the command rate r^: 



where 

(e n /A, e^/A) , if (a^ is chosen 

(e I3 /A, e J4 /A) , if (bj) is chosen , 

(e m /A, e n2 /A) , if (a^) is chosen 

( r 3 * r i) " { 

(e n3 /A, e E4 /A) , if (b n ) is ohosen , 
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r ' max * "“(‘'i 1 ' |r 2' • |r 3' > |r 4 l) ’ 


(17) 


and r is the maximum limit of the gimbal angle rate, 
g 

Remark 4 ; r.' is the gimbal rate of the i-th CMG required to realize the 

desirable gimbal angles obtained in Step 3. The meaning of equation (16) is that 
if any of the r'. T s exceed the maximum limit of the gimbal rate r (which is 

determined by CMG hardware) , then gimbal rates r 1 . ~ r^ are limited propor- 
tionally to minimize the effect of undesirable output torque due to the hardware 
restriction. 

This completes one cycle of calculation of the command rate r^ from 
a (n) and t^ (n) . At time (n + 1) A, upon obtaining the new state a (n + l) , a 
new cycle begins from Step 1. 


MOMENTUM DISTRIBUTION SCHEME 

In this section a momentum distribution scheme between h Il andh m’ 
i, e. , a specification of a pair of functions f^. and f^, is given. 

Let x and x be 
1 £ 


X 






(18) 


These values express the capability of the two pairs in producing momentum in 
the direction of the Y-axis when h^ = h^, ^112 = ^d2" these variables, 

we specify that 
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tjlh,. <*(»)> = Vd3/<V X 2 ) + g 


(19a) 


and 


‘II 


th^ a (n)] = 


X 2 h d3 /(X 1 +X 2 ) ' g 


(19b) 


where the first terms in the right-hand sides denote a proportional distribution 

of h , and g of the second term represents an additional (or excessive) distri- 
d3 

bution. Notice that for any g, f and f II satisfy 


f I + f n h d3 


Therefore, we can select any value for g. 1 


First, a desirable g, g* will be given as a function of h^, an d x^. 

Then the value of g will be given as a function of g* and the present state a (n) . 
This two- stage approach is taken to make the present additional distribution g(n) 
converge to the desirable additional distribution g* without too much change for 
each sampling interval. In the following, g* is given. First let us define two 
candidates for g + : 


g a = 


X 1 X 2 


'90 h 


d3 


'90 h 


d3 


0.9cos \x l + x 2 J+ ('n/T- 0.9) cos + x g 


(20a) 


and 


X 1 X 2 


'90 h 


0. 8 cos 


d3 


X 1 + X 2> 


(20b) 
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We then select one of g and g as g* in the following way: 

£L 0 


g 


* = 


where 


x 2 x 

v “ g c - rr^r h d3 8 h d3 8 e c + rrir h d3 
12 12 


g^, otherwise 


( 21 ) 


(0. 5 + k ) g + ( 0. 5 - k ) g , if g* •= g at the last cycle 

13- It) 3 


g c = « 


( 22 ) 


( 0. 5 - k ) g + ( 0 . 5 + k ) g , if g* = g, at the last cycle 
1 3 lb b 


and 0 g k^ 2 0. 5 (for example, k^ = 0.2). 

A desirable momentum distribution | f^* , J which corresponds to g* 
is given by 


V = x iV l V , i ) + 6,1 


(23a) 


and 


f * 
II 


X 2 h d3 /(X 1 + 


X 2 ) " S* 


(23b) 


Figure 5 shows a schematic diagram of f^* and f^* as a function of h for a 

given pair of x and x . 

12 


♦ 
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a. If g* = g at the last cycle. b. If g* = g, at the last cycle. 

cl D 

Figure 5. Schematic diagram of the desirable momentum distribution. 

Now let us specify a relation between g and g*. The additional distri- 
bution g(n) at the present state a (n) is given by 


g(») = h n ( n ) 


x x (n) th n (n) - h II1 ( n) ] 
X 1 + X 2 ^ 


(24) 


where 


x ! ( n ) £ 


- \i2 2(n) 


and 


16 


<«) 




4 h , 


- h 


12 


(n) 


Using g(n) , we specify g as follows: 


[g» , u | g* - g( n) I S g max 


e = 


(25) 


g(n) + g max sgn(g*-g(n)l , If Ig* - g(n) I > S max 


where 


g 


max 


k A r h ir/180 
2 g * 


(26) 


and k is a constant (for example, k = 0. 5) . 

A ** 

The above equation means that, although the additional distribution g at 
the next sampling instant should be as close to g* as possible, its change in one 
sampling interval should not be larger than a maximum allowable change g max * 

The constant k determines a degree of instability of the singularities. There- 
A 

fore, if the torque command is kept zero, any initial state converges to a 
desirable momentum distribution state in a finite time. 

Returning to discussions of the desirable momentum distribution as was 
shown in Reference 2, there are two types of internal singularities (i.e. , r .i' / 
singularities inside the momentum envelope) : 


(a) - - a 2 = ±90 deg or = - <* 4 = ± 90 deg 


and 



oi 2 1 = l a 3 ■ a 4 l = 180 deg . 


(b) !a 1 - 


When the OMEGA steering law is used, there are stable singular states 
of type (a) although those of type (b) are unstable. 

Singular states (a) and (b) are represented by lines a , a , and point b 

X 

in Figure 6. From the previous discussions it is clear that, as the time evolves, 

points on lines a and a (including point b) are transferred to points on the lines 
X & 

of f * and fjj* (for example, if the torque command is zero, point b will be trans- 
ferred to point b 1 in the figure) . Therefore, all the internal singular states are 
unstable when the proposed steering law is used. This instability causes quick 
movement out from singular states and will make very small the possibility of 
facing unfavorable torque command just at the singular state. Moreover, even 
if an unfavorable torque command is applied just when the system is in a sin- 
gular state, the control deterioration is very small. This will be shown by 
digital simulation results in the section, Computer Simulations. 

A realistic situation in which serious control deterioration may occur will 
be the passage of the system state through a neighborhood of a singularity caused 
by a fluctuating torque command. This fluctuation is unavoidable because of 
various disturbance torques and sensor noises. A model of such a situation is 
shown in Figure 7, where a small sinusoidal torque command in the YZ-plane 
with a little bias in the Y-axis is applied to the system at zero momentum state. 
This sinusoidal torque command forces the system to pass near, or hit, a 
singularity. 

Now let us show by an intuitive argument that the hysteresis introduced 
at two jumps of f * and f * keeps the control deterioration very small in this 

situation. In Figure 5 the passage of the system close to a singular state cor- 
responds to the passage of the value of h close to one of two jump points. 

Uu 

According to the fluctuation of torque command, the values of x and x wili 

1 2 

also fluctuate. If there is no hysteresis [that is, k^= 0 in equation (22)] , this 
fluctuation may cause a fluctuation of g* between g and g , which causes control 

cl D 

deterioration caused by hardware limit on the gimbal rate. By the introduction 

of hysteresis, however, most of the fluctuation of g* between g and g could be 

a b 

avoided, thus keeping the control deterioration to a minimum. 
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Figure 6. Relation between internal singular states and the 
desirable momentum distribution. 

Because of these two features of this steering law, the instability of the 
singular states and the hysteresis of the desirable momentum distribution, we 
can say the following: If any control deterioration occurs because of a sin- 
gularity, keep the torque command zero, and in a short time the trouble will 
be gone. 

Candidates for g*, g and g, of equation (20) have been selected as a 

a d 

simple pair which satisfies the following requirements: 
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0 


TIME 


cx 





TIME 


a. Torque command. 


X 



Figure 7. A model of fluctuating torque command and 
corresponding transition of total momentum. 

1. The value of g should be \f2 for x = x = 2 and h = 0 in order to 
a 1 2 d3 

give the gimbal angles {a. = a = 45 deg, a=a =-45 deg} for the zero momen- 

turn stationary state. These gimbal angles are desirable for the LST because 
the possible output torques in the Y and Z directions are balanced. (A more 
detailed discussion is given in the appendix. ) 
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(27) 


2. 0 = £ b = = 


min (x , x 2 ) 


TT < x i + V lh d3 l) 


1 2 


for 0 ^ x i < 2 , 0 i x 2 i 2 , Ih^l = X x + x 2 , x^^ + x g * 0 . 


3. g and g, should be smooth. 

°a b 

So far we have neglected the case of | h | 5 2h*, Ih^l = 2 h *> or 
x (n) + x (n) = 0, which may cause trouble in treating equations (19) through 

1 dt 

(21) and (24). This trouble can be easily avoided, for example, by using, 
instead of equation ( 1 8) , 


J 4h 2 - h.„ 2 , if 4h “ - h . “ £ e 


* dl 


*1 = 


1 ’ 


2 2 2 
* - h di = e i 


2 2 2 
if 4h, - hfc < £l 


"^ 4h * 2 - h dl 2 > 11 4h *" - h ^‘ * £ 


*2 - 


£ 1 ’ 


2 2 2 
* - h d2 = 


2 2 2 
if 4h. - h d2 < « 


where e is a small positive number (for example, = 0.0001) and by using, 
instead of equation ( 25) , 

x i (n)[h i2 (n) - h n2 (n)1 

g(n) = h f2 n - Xl (n) + x 2 (n) + £2 

where e is a small positive number (for example, e = 0*00001) . 

2 2 
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CMG-OUT OPERATION 


We assume that a signal is sent to a control computer continuously from 
each CMG indicating whether it is operating. We will consider the case with 
only one CMG-out. Under this assumption, we will develop a steering law which 
automatically transfers the system to the state of CMG-out operation when a 
CMG fails and, moreover, automatically resumes the state of normal operation 
when the out-CMG recovers. 

When a CMG is out, there is no capability for an arbitrary momentum 
distribution. For example, when the first CMG failed, h =_h Q . Hence, the 

*“l Li 

magnitude of h is fixed to h . 

— l * 

Since h T „ should take the value h , h T „ can take only one of the two 
12 dl II 

T 2 2~ 

values ± \I h^ - h ,, . Furthermore, the direction of h T should not be far 
* dl —I 

away from that of the present momentum h^ ( n) . A simple way to achieve this 

may be given as follows: 

’ 

+ ^ h * 2 " h dl 2 ’ if h Ii (n) > 0 

h n = 

or if h^(n) = 0 and 

- '^7 ~ h dl 2 ’ if h n (n ) < 0 

or if h ^ ( n) = 0 and 


d3 


d3 


< 0 


h m " ' h d 3 + h n- 
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Taking into cons ids rat ion that the total duration of the CMG— out opera 
tion is (and should be) much less than that of normal operation and that the 
operation should go back to normal when the out-CMG recovers, the CMG- out 
operation should be performed by a steering law which is close to that of 
normal operation, even if that law has some redundancy. 

From this point of view, the following modification to the steering law 
for normal operation seems to be good to include the CMG-out operation. 

Let w = [w ! w ! w j wJ T be a CMG-out signal vector. When the i-th 

CMG is not out, w t = 1 and when it is out, w. = 0. 

In step 1, instead of equation ( 6) , we adopt. 


h (n) - [w x cos « 1 (n) + w g cos « 2 (n)] , 


h^ (n) = h* 1^ sin (n) + w 2 sina 2 (n)] , 


h m ^ = h * 1 w 3 COS a 3 + w 4 COS a 4 ^ ’ 


and 


h H2 ^ = h * [w 3 sino! 3 + W 4 Sina 4 ^ * 


In step 2, concerning equations (7) and ( 8) , three cases should be 
treated separately: 

Case 1. w + w 2 + w^ + w^ = 4 
( normal operation) . 
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1 


Case 2. w + w = 

1 w 


(a CMG of the pair I is failed) . 


Case 3. w + w = 1 

O 4 


(a CMG of the pair II is failed) . 

In Case' 1, use equations (7) and (8). In Case 2, use 


h Il 


{sgn[h I:L (n)l + [ 1 - sgn [h^ (n)] sgn (h d3 )} '/h* 2 - h dl 2 . 


0, 


lfh dl < 1 


“ h dl * 1 


and 


h m = - h d 3 + h u • 


In Case 3, use 


h m = 


{sgn [hm (n)l + tl - sgn [h ni (n)] sgn (h, Q )} 


/. 2 2 
d3 /J ^ h * h d2 ’ 


0, 


if h < 1 
d2 


ifh d2 8 1 


and 


h n h d3 + h m * 
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Steps 3 and 4 need no modification. In step 5, instead of equation (17) , 
we adopt 


. r' 


max 


•(Wjiqi, w a lrj!, Wglr^l, w 4 |r'J) 


COMPUTER SIMULATIONS 


A program for computer simulation of the CMG system with the proposed 
steering law was made using BASIC language and a minicomputer. Values of 
constants were selected as follows: 


j3 .= 30 (deg) , h* = 1 (normalized) , 

A = 2 (sec) , r = 2 (deg/sec) , 

S 

k = 0.2 , k = 0. 5 , 

1 ^ 

e = 0.0001 , e = 0.00001 . 

1 4 


Some results of the computer simulations are given in Figures 8 through 
12. Figure 8 shows the transition of the gimbal angles o', and the output torque 

in the XYZ-coordinate from the zero momentum stationary state, a = 

[45 deg j — 45 deg [ 45 deg j — 4 5 deg] , when a torque command of magnitude 
0. 01 (normalized number by h^ ) is applied in the direction of the Y-axis. No 
loss or deterioration of control occurs passing through the singularity as in the 
case of the OMEGA steering law. Moreover, this singular state is unstable, as 
shown in Figure 9, where the torque command was changed to try to stop the 
system at the singular state but the attempt failed. 

The effect of unfavorable torque commands at a singularity is shown in 
Figures 10 and 11. Figure 10 shows the case where the torque command of 
magnitude 0. 01 in the direction of the Y-axis was suddenly changed to the 
direction of the minus X axis ( most unfavorable direction) just when the 
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Figure 8. Transition of girabal angles and output torque 
for Y-axis torque command, 
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Figure 10. Effect of unfavorable torque command at singularity 
(command magnitude - 0,01) , 
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Figure 11. Effect of unfavorable torque command at singularity 
(command magnitude = 0. 005 ) . 
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Figure 12, Response to a fluctuating torque command. 


30 


system reached the singularity. The figure shows a delay of 12 sec in response 
of the output torque to the command torque. Figure 10 shows the case where the 
magnitude of unfavorable torque command is 0. 005, i. e. , one-half of that in the 
previous case. The delay in this case is 6 sec. Generally, the smaller the 
magnitude of unfavorable torque command, the smaller the delay in output torque. 

Figure 12 shows the response to a fluctuating torque command which is 
composed by adding an X-directional sinusoidal component to that in Figure 7: 


t 

— c 


(t) = 


0. 0025 sin (10 r) 

0.005 sin (5 t) + 0.0005 
0. 005 cos ( 5 t ) 


T 

The initial state is a = [21.36 deg I - 21. 36 deg j 81. 88 deg ! - 81. 88 deg] , 

i 1 l T 

which corresponds to h^ = [0 J 1. 52 J 0] . 

If no hysteresis is introduced into the desirable momentum distribution, 
that is, if k = 0 in equation (22), the response to the same torque command 

becomes that given in Figure 13. It is clear by comparing Figures 12 and 13 
that the introduction of hysteresis keeps the period of control deterioration to a 
minimum (it should be noticed that slewing of and for 180 deg cannot be 

avoided) . Any steering law for the roof-type configuration with no hysteresis 
will suffer to a certain extent from such a control deterioration as that in Fig- 
ure 13. Figure 14 is a typical example of operations with a CMG failed where a 
torque command of 0.01 is applied in the direction of the Z-axis. As seen from 
the figures, the control performance is quite satisfactory. 


DISCUSSION 


Because of Step 4, our steering law redistributes away from undesirable 
initial states in a very fast way as shown in Figure 15, where the initial gimbal 
angles are (0 deg, 180 deg, 90 deg, - 90 deg) in (a) and (- 30 deg, - 20 deg, 
-135 deg, 45 deg) in (b). Although slightly undesirable output torques are seen 
in the figure, they are negligible. 
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Figure 13. Response to a fluctuating torque command 
(case of no hysteresis in desirable momentum distribution). 





Even if the length of the sampling interval is very large, the control 
performance of our steering law has no theoretical error as far as the control 
purpose stated in the section, Roof-Type CMG System, is concerned, except 
when it is impossible to meet the torque command because of the gimbal rate 
limitation or momentum saturation. This might be an advantage if changes of 
sampling interval during the operation of the CMG system are desired for some 
reason. It should be noted that a longer sampling interval may cause a control 
deterioration due to fewer updatings of the torque command. 

In most steering laws developed so far, additional logic is necessary -to 

prevent the gimbal angles from getting into an oscillation with maximum rate 
after reaching any momentum saturation. When our steering law is used, there 
is no such oscillation. Roughly speaking, the system moves toward a saturation 
state where the direction of the total momentum is the same as that of the 
desired total momentum. We can also make the system stop whenever it reaches 
a saturation state, and keep that state until a torque command in a direction of 
desaturation is applied, simply by adding the following additional step to our 
steering law. 

Step 6. 


r - 0, i - 1, 2. 3, 4 

ci 


if 


h I > < W l + W 2> h * 


(« r -«2) Hr 1 tw 1+ |r 2 lw 2 + w 1 w 2 ] = 0 , 

; ■ ■ ■ k 

! are satisfied, or if - 
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h n > (w i + w 2 } h * 

and 

(a 3 - « 4 ) C I r 3 1 w 3 + k 4 tw 4 + w ? w 4 l = 0 , 


are satisfied. 


So far, we have discussed our steering law with one specified desirable 
momentum distribution {f^, f^*} given by equation (23) . The answer to the 

question of what the desirable momentum distribution is might vary, depending 
on each engineer and on the capability of the CMG hardware used. One of the 
features of our steering law is that, since our momentum distribution scheme 
treats the momentum directly, it can easily realize any reasonable momentum 
distribution. To illustrate this, an example is given in the following. In this 
example, we attempt to achieve a control performance similar to that of the 
OMEGA steering law. The desirable additional distribution g* is selected as 


X 1 X 2 

g* = — cos 


90h d3 
X l + x 2 


A schematic diagram of the corresponding desirable momentum distribution is 
given in Figure 16. A simulation result for the case of unfavorable torque 
command at a singularity is given in Figure 17. The figure shows a delay of 
44 sec. Several other simulation results have also shown a control perform- 
ance which is quite similar to that of the OMEGA steering law. 

The computing requirements will now be discussed. According to the 
author’s experience, the length of the program for our steering law is about 
1.2 times that for the OMEGA steering law when a proportional limit on rate 
command [like equations (16) and (17)] is included in the OMEGA steering 
law. Since our steering law is given in five steps, each of which has a very 
simple physical meaning, it will be easy to modify or to simplify it in accord- 
ance with various requirements in practical applications without changing its 
main concept. When this point is considered, the computing requirement of 
our steering law is very reasonable. 
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It will be possible to give a control performance similar to that of our 
steering law by adding the same idea as that in the section, Momentum Distri— 


bution Scheme, to the OMEGA steering 
law. However, the author believes that 
the approach taken in this report is 
easier to understand and has various 
advantages over any modification of the 
OMEGA steering law. 


CONCLUSIONS 


A steering law is proposed for 
a roof-type configuration of the 
SGCMG system which is obtained 
by regarding the CMG system as a 
sampled-data system and providing a 
new momentum distribution scheme. 
This scheme is designed to bring any 
state of the system to a state with a 
predetermined desirable momentum 
distribution, which has two jumps with 
hysteresis around singular states. It 
is analytically shown that these jumps 
make all the singular states unstable 
and that these hysteresis effects make 
the system relatively insensitive to 
singularities. 



Figure 16. Schematic diagram of desir- 
able momentum distribution correspond- 

/(v 


ing to g* = x x 2 cos [ 90 h^ 


* 2 )l/4. 


With these two features it is 

expected that the steering law will give a control performance which is good 
enough for practical applications. Results of the preliminary computer simu- 
lations entirely support this expectation. 
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APPENDIX 


DES I RABLE ZERO MOMENTUM G IMBAL ANGLES 
AND TORQUE ENVELOPE 


The shape of the total momentum envelope has been the criterion for 
determination of the skew angle j3. However, when a CMG system is accom- 
panied by a device for desaturation such as a magnetic torquer in the case of the 
LST, most of its operation period would be spent around the zero momentum 
stationary state. In such a case, the shape of the envelope of the possible torque 
output at the zero momentum stationary state will also be important. If this 
shape is not proper, it may happen that a much larger rate command is necessary 
to give an output torque with a fixed magnitude in one direction than that necessary 
in another direction. 

This torque envelope at the zero momentum stationary state could also 
serve as a criterion for determining desirable zero momentum gimbal angles, 
whereas the momentum envelope is not a proper criterion for this. 

In this appendix it is shown that, for the LST, 0 = 30 deg, which was 
recommended from the viewpoint of the momentum envelope, is reasonable also 
from the viewpoint of the torque envelope, and that desirable zero momentum 
gimbal angles are = 45 deg, = - 45 deg, = 45 deg, a ^ = - 45 deg. 

First, all sets of gimbal angles which give zero momentum will be 
obtained. From equation ( 3) , 


h I2 ' ° ' 


h II2 = ° • 


and 


h Il " h IIl 
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Hence, any set of gimbal angles for zero momentum should satisfy 


a 


1 



( A-l) 


or 


l a i _Qt 2* = ' a 3 _a 4' = 180 de S ♦ 


(A-2) 


I i T 

Second, a torque envelope will be defined. An output torque [t ( t 1 1 ] 

' x , y | z 

at any state a and gimbal rate & is given by 


t = d h /dr = sin B (a, cos a, + a n cos a n + cos a_ + a . cos a .), 
x x 112 2oo4 4 


t = d h /dT = - (o' sin a. + a n sin b 0 + a 0 cos a 0 + a . cos a ) , 

X V 1 1 i £ O d 4 4 


and 


t = dh /dr = cos j3 (a. cos a, + a 0 cos a„ - a 0 cos a - a. cos aj. 

(A-3) 

A torque envelope at state a is defined as the envelope of all possible output 
torque vectors under the constraint |a^| = 1. 

Third, the torque envelope for zero momentum states will be obtained. 
For a set of gimbal angles satisfying equation (A-2) , the torque envelope Is a 
diamond- shaped plane (Fig. A-l) and this state cannot produce an output torque 
perpendicular to this plane. This corresponds to the fact that this state is 
singular. Hence, equation (A-2) is not desirable as a zero momentum state. 

For a set of gimbal angles satisfying equation (A-l) , the torque envelope is a 
duodecahedron shown in Figure A-2 where 
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4 sin /3 | cos 1 , 

(A-4) 

4 ! sin a \ , 

(A- 5) 

4 cos (3 | cos a | ; 

(A-6) 


t *, t * and t * are interpreted as the 
x y z 

normalized maximum possible output 
torques in the X, Y, and Z directions. 

It can be assumed that 0 5 = 90 deg 

without any loss of generality because 
of the symmetry of the system. 

Now, the question of what values 
should be selected for (3 and to make 



the torque envelope have a good shape 

will be discussed. For the LST, the 

principal moments of inertia 1,1. 

x y 

and I satisfy I si £61 . It is con- 
z y z x 

sidered to be reasonable to make t *, 

x 

t * and t * proportional to I , I and I . 
y z x y z 


at 

Figure A-l. Torque envelope at state 
a = (30 deg, -150 deg, 30 deg, 

-150 deg). 


In order that t * = t^* should be satisfied, from equations (A- 5) and 
(A-6), and /3 should satisfy 


ot 


1 


tan 


-1 


(cos (3) 


(A-7) 


Then, 


t * 

y 


t * = 4 sin [ tan 
z 


(cos /3) ] 


(A- 8) 
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0 






Figure A- 2. Torque envelope at the state 
^l = -° i 2 :=±0! 3 = :F “4 ) - 

On the other hand, from equations ( A-4) and ( A- 6) , the following relation 
is obtained: 

t x * = t * tan/3 ; (A- 9) 
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a t *, t * and t *, given by equations (A- 7) , A-8), and(A-9), are plotted 
1> y z x 

in Figure A-3 as functions of the skew angle /3. 

The value of B for which t * = t * = 6 t * is satisfied is B ~ 9« 7 deg. But 

y z x 

the decrease of t * and t * due to the increase of j3 from 9. 7 deg to 30 deg is 

y z 

only 6. 7 percent, while the increase of t * is more than 200 percent. Moreover, 

the angle 30 deg is geometrically simple. Therefore, the best skew angle would 
be 30 deg from the viewpoint of the torque envelope. For j3 = 30 deg, the angle 
given by equation (A-7) is 40. 893 deg and 


t * = t * = 2.619 , 

y z 

( A-10) 

t * = 1.512 . 

x 


The values of t *, 

y 


t * and t * for fl = 30 deg and a = 45 deg are given by 
z x 1 


t * » 2.828 , 

y 


t * = 2.450 , 

z 


t * e 1.414 
x 


These values are not too different from equation (A- 10) . The angle 45 deg is 
also a geometrically simple angle. Hence, = 45 deg, = - 45 deg, = 45 

deg, and = - 45 deg are recommended as the best zero momentum stationary 

gimbal angles. 

It can also be said from the veiwpoint of the torque envelope that an 
angle much larger than 45 deg is not desirable as For example* t *, 

t *, and t * for B = 30 deg, a, = a = 60 deg,' a. = a = -60 deg are given by 
Z X 13 2 4 
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t * = 3.464 , 

y 

t * = 1. 732 , 

z 


t * = 1.000 . 

x 


Roughly speaking, this means that, for a torque command in the Z-axis, the 
gimbal angles should be driven with a speed twice that for a torque command in 
the Y-axis with the same magnitude. This, of course, is not desirable. 
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